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Abstract 

The chiral symmetry is applied to an extended constituent quark potential model. 
With random phase approximation (RPA), a small component effect is added to 
the constituent quark model. To obtain the pseudoscalar tt meson as a Goldstone 
boson, the quark effective potentials are modified in the model to account for the 
dynamical breaking of chiral symmetry. Also the vector p meson is calculated and 
the KSRF relation about tt and p meson decay constants is derived in the model. 
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1 Introduction 

Due to the complication of a non-Abelian SU(3) gauge theory, quantum chro- 
modynamics (QCD), which describes the strong interaction, has many non- 
perturbative features such as the dynamical breaking of chiral symmetry and 
quark confinement. In the study of hadron structure at low energy scale, two 
kinds of models of QCD are often used, each incorporating some important 
QCD features. The constituent quark potential model which incorporates the 
QCD quark confinement has been impressively successful in hadron spec- 
troscopy and decays [1], except for the pseudo-scalar tt meson which is a 
Goldstone boson and has very low mass. It was shown recently that the RPA 
y-component must be considered in the treatment of pions [2]. On the other 
hand, the Nambu-Jona-Lasinio (NJL) model [3] and several of its extensions 
[4,5,6] which incorporate the chiral symmetry can describe the it meson very 

Email address: dwz@pku.edu.cn (Weizhen DENG). 



Preprint submitted to Elsevier Science 



1 February 2008 



well as a massless Goldstone boson of dynamical breaking of chiral symme- 
try. However, the NJL model lacks the important QCD property of quark 
confinement. 

In Ref. [7], we proposed an extension to the constituent quark model to com- 
prise the RPA ^/-component quark excitations. As in NJL model, we start from 
an effective quark Hamiltonian. For meson structure, this leads to a Hamil- 
tonian with two channels — the ordinary valence quark x-channel and a new 
^-channel, which is an extension of the quark potential model. This extended 
quark potential model has a new coupling potential connecting the ordinary 
valence quark a:-channel and the small component ^/-channel. The coupling 
potential can be related to the ordinary valence quark potential in the model. 

In this paper, we implement the chiral symmetry in the extended constituent 
quark potential model. First we introduce the extended quark potential model. 
In Sec. 3, we analyze ir and p mesons using this framework. The restriction 
from chiral symmetry is considered. With the meson wave functions containing 
^/-channel quark excitations, KSRF relation [8,9] about 7r and p meson decay 
constants is derived. Finally, we summarize our discussion. 



2 Extending Quark Potential Model with RPA 



QCD perturbation theory is very successful when applied in the high energy 
processes. In the low energy regime, Hamiltonian approach can provide an al- 
ternative method to understand the non-perturbative features of QCD, such as 
the quark confinement. Many efforts were made to obtain an effective Hamil- 
tonian from the exact QCD Lagrangian [10,11,12]. Here we start from the 
effective quark hamiltonian, assuming the quark fields have been separated 
approximately from the gluon fields 



H = H 2 + H 4 , (1) 
H 2 = J d 3 x^(x)(-ia • V + Pmo)V(x), (2) 

#4 = ^E / d 3 xd 3 yK e (x - y)*(x)rM*)*(v)rMv), (3) 

where H 2 is the free current quark Hamiltonian and itlq is the current quark 
mass, if 4 includes all effective quark interactions. K a {x) is the kernel function 
for coupling vertex T a . 

Apart from a constant, the Hamiltonian (1) can be written in normal order as 
H = H +:H A : (4) 
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where H represents the free particle energy of constituent quark and its gen- 
eral form is 



H 



J d 3 x^\ 



X 



-L4(-V 2 )tx-V + pB(-V 2 ) V(x), 



(5) 



where the dependence on the three-vector momentum arises from the instan- 
taneous approximation of quark interactions. 

The two functions A{k 2 ) and B(k 2 ) should be obtained from a self-consistent 
mean field calculation. However in a constituent quark potential model, the 
quark interactions contain the linear quark confinement which has severe in- 
frared divergence in momentum space. On the other hand, the constituent 
potential quark models always assume that the constituent quarks have fixed 
masses. As an extension of the constituent potential quark model, here we 
assume that the A and B functions can be approximated as constants. So we 
treat the constituent quark mass as a parameter, m c = B (with A = 1). The 
free particle Hamiltonian becomes 

H = [ d 3 x^(x) l-ia • V + (3m c \ (6) 



In this paper, we will only consider meson structures. In naive quark model, 
meson fields are approximately expressed in terms of local quark fields as 

4>m{x) = a M *(x)r M *(x), (7) 



where Ym = ilb^i for iso- vector pseudo-scalar it mesons or Ym — 7 • tiTj for 
iso- vector vector p mesons. According to 1/N expansion, the above meson 
structures are exact in the limit N c — > oo [13,14]. 

In the non-relativistic limit, as quark momentum k <^ k°, one has [15] 



<p M {.x) = a M ^2 J 



d 3 h d 3 k 2 1 1 
(2^)3 J2^2kj2kl' 



-i(ki+k 2 )-x 



S1S2 

x [u(-fci, s 1 )T M v(-k 2 , s 2 )b' t (-k l , s 1 )d' t (-k 2 , s 2 ) 
+ v(k x , si)T M u(k 2 , s 2 )d(ki, s 1 )b(k 2 , s 2 )] , 



(8) 



for 7T and p mesons. On the other hand, the field operators of n and p can also 
be expressed in the Q and Q} of meson creation and annihilation operators 



r d 3 k 1 

M*) = J 7^)3^0 [Q M (k)e-^ + QUk)e^ 



(2tt) 3 2k° 



(9) 
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If the mesons are composed of creation of valence quark pairs as in the poten- 
tial models, one deduces from the equivalence of eqs. (8) and (9) 



xw(fci, s 1 )T M v(k 2 , s 2 )b t (k 1 , s 1 )d\k 2 , s 2 ). (10) 
However, generally one can only obtain 

x [X M iZ(fei, si)r M w(fc 2 , s 2 )6 t (fci, Si)d f (fc 2 , s 2 ) 
+ YM^(-fci, si)r M M(-fc 2 , s 2 )d(-fei, si)6(-fc 2 , s 2 )\ ■ (11) 

This is just the RPA excitation operator — the meson state is a superposition 
of creation and annihilation of qq pairs on the vacuum. The second part, i.e. 
the y-component is simply discarded in the ordinary quark potential model as 
a small component. 

Thus, as an extension to the potential model, we take the mesons as excitation 
modes of the vacuum of the RPA type: 

\Q) = Q ] |o> . (12) 



For a meson in the rest frame, after we take account of the quark interactions, 
excitation operators are 



r d 3 k 1 

^ t = Ey (2^)3 2fc° [ ;r ( fc ' s i' s 2)& t (fcs 1 )rf t (-fcs 2 ) 



+y(k, s 1 ,s 2 )d(-ks l )b(ks 2 )\ , (13) 
where the time reversals are defined as 



b(ks) 
d{ks) 



-iy/ 2 - s b(k-s), 

-l) 1/2 - s d(k-s), 



(14) 
(15) 



and x and y are the RPA amplitudes determined from the well-known RPA 
equation of motion[16] 



(0| [6Q, [H, Q 1 ]] \Q) = E(Q\[6Q,Q^]\Q). 



(16) 
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We obtain the equations for RPA amplitudes 



(E-2k°) x(k, Sl ,s 2 ) 
f d 3 k' 

: 2 / 7^3 s u s 2 \U\ k', s 3 , s A }x(k', s 3 , s 4 ) 

+(k, si, s 2 | V | fc', s 3 , s A )y(k', s 3 , s 4 )] , (17) 
(s + 2A; )y(fe, Sl , S2 ) 



= - S / K fc > s i> s 2 I V | fc', s 3 , s 4 )x{k', s 3 , s 4 ) 



+ (fe, si, s 2 I ^ I fc', «3, s 4 )y(k', s 3 , s 4 )] . (18) 

Here we have introduced two effective potentials [/ and V. Their matrix ele- 
ments are 



1 1 

(k, si, s 2 | C/ | fc', s 3 , s 4 > = -7^7^ E " 

M(fcs 1 )r e M(fc's 3 )t;(-fc / S4)r e t;(-fcs 2 ), (19) 

(fe, si, s 2 1 y | fe', s 3 , s 4 ) = E - 

M(fcsi)r e fJ(-fc / s 3 )S(fc / s 4 )r e t;(-fcs2). (20) 



Thus, with the RPA approximation, the quark potential model is extended to a 
coupling system with two channels. The first part, i.e. x-amplitude x(k, si, s 2 ) 
which we will call the x-channel wave function, is the large component which is 
the sole valence quark contribution in the ordinary quark potential model. The 
second part, i.e. the y-amplitude y(k, Si, s 2 ) which we will call the y-channel 
wave function, is the small component which is discarded in the ordinary 
quark potential model. The potential U interacts only within each individual 
channel while the new coupling potential V couples the two channels together. 
If the coupling between the two channels V is small, the system decouples. 
The extended potential model reduces to the ordinary potential model and 
the potential U is just the ordinary quark potential. Here we will show that 
the coupling potential V is important due to chiral symmetry . 

One can use the Dirac kets and bras to simplify notations. x(fc,si,s 2 ) and 
y(k, Si, s 2 ) are represented by two kets \x) and \y) respectively, 



x(k,s 1 ,s 2 ) = (k,s 1 ,s 2 \ x), (21) 
y(k,s 1 ,s 2 ) = (k,s 1 ,s 2 | y). (22) 

Eqs. (17) and (18) can be written concisely as 
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(P-2A; ) \x) = U\x) + V\y), 
(E + 2k°) \ y ) = -V\x)-U\y), 

or the familiar matrix form 



(23) 
(24) 



2k° + U V 




\x) 




\x) 






= E 


-V -2k° - U 




\y) 




\y) 



(25) 



where k° = \jm 2 c + k 2 . 



By definitions (19) and (20), the potentials U and V will be evaluated from 
the same kernels in the quark interaction if 4, thus they can be related to each 
other. In the non-relativistic limit k, k' — > 0, we have 



U — S — V — A(Ti ■ <T 2 — T<7\ • cr 2 , 
V = P- Vct x -o-2-A- Tct x ■ <7 2 , 



(26) 
(27) 



where on the right hand side of equations, the kernels are divided into five 
types as in ref. [17], i.e., scalar (S), pseudo-scalar (P), vector (V), axial- vector 
(A), and tensor (T); cr's are the pauli spin matrices. 

To normalize the model wave functions, we apply the standard boson commu- 
tation relation 



IQ(P), Q\P')\ = (27r) 3 2P°5 3 (P - P') 



(28) 



to rest mesons. We obtain 
(Q\Q) = (2tt) 3 2P5 3 (0). 



(29) 



Then we find the normalization relation of RPA type 



Sl,S2 ' 



(30) 



Again with Dirac kets and bras, its form is 
2E=(x\x)-(y\y). 



(31) 



With the normalized meson wave functions, one can calculate the static prop- 
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erties of mesons. The tt weak decay constant is denned as 

iUWa = (0 | ^(0)7^75^^(0) | n 3 (q)). (32) 



For pions at rest, the above equation becomes 

1 



ifirmJtj = (0 | *(0) 7 75-r i *(0) | Q nj ). (33) 



i 



In the non-relativistic limit, we obtain 

iUm* = (r = 0, S = 0, M s = | x) + (r = 0, S = 0, M s = | y), (34) 

where (r, S = 0, Ms — | x) and (r, S* = 0, Ms = | are the coordinate 
wave functions in x- and y- channel of the n meson with its total spin S = 
0, M s = 0. The p electro-magnetic decay constant is defined as 

m? - 1 

= (0 I *{0)-f-TM0) I PM), (35) 

where i, j are the isospin indices, and k is the p meson polarization index. For 
static p mesons 

77l 2 - 1 

ytfy = (0 I *(0) 7 fc -T^(0) I Q pjk ). (36) 

In the non-relativistic limit, one obtains 

2 

in 

_£ = ( r = 0,5 = 1,M 5 I x) + (r = 0,S = 1,M S | y). (37) 



3 7r and p meson Properties 

Following the quark potential model [1], we choose the effective quark inter- 
action H A as a vector interaction plus a scalar confinement interaction 

Ha = \S d3xd3 v K s( x ~ y)y(x)v(xMy)*(y) 

+ l - J d 3 xd 3 yK v (x - yM*b M y * (aW^y *(y)- (38) 
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According to lattice calculation, the scalar confinement kernel is linear in r 



K s (r) = br. (39) 



The vector kernel is taken from the one-gluon exchange plus a constant 

K v (r) = c+^1, (40) 



where a s (r) is the QCD running coupling constant. 

In our analysis, we use the simple harmonic oscillator wave functions for pions 
and rho mesons as a first approximation. In this way the underlying physics 
is emphasized. 



x(k,s 1 ,s 2 )=X* 000 (k){^s 1 ^s 2 I SM S ), (41) 
y(k,s 1 ,s 2 ) = Y* 000 (k)(±s 1 ±s 2 | SM S ), (42) 

where S, Ms are the spin quantum numbers of the meson, and ^oooW is the 
ground state wave function of harmonic oscillator in momentum space 



In the non-relativistic limit, the two potentials U and V are 

U = br+ (c + ^Fx-Fa, (44) 
V = ctr 1 ■cr 2 F 1 • F 2 , (45) 

where 



Fi = ( 



for quarks, 

(46) 



A c A* 

for antiquarks. 



2 2 



Here, one must be careful with the one-gluon exchange interaction — . It arises 
from the covariant form in momentum space, where Q = k — k' is the 
4- momentum of the exchanged gluon. For potential U, k — (k°, k) and k' = 
(k'°, k') are the 4-momenta of incoming and outgoing quarks respectively. Thus 
in non-relativistic limit, l£ ~ ^ which is a Coulomb potential in coordinate 
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space. However for coupling potential V, k — (k°, k) and —k' = (k'°, —k') 
are the 4-momenta of qq quark pair which were created or annihilated on 
the vacuum. So in non-relativistic limit, ~ ^ ~ 0. The term — is 
suppressed in V. 



It is well known in NJL model, the chiral symmetry can be broken dynami- 
cally with a nonzero vacuum condensate coming from quark interaction. The 
current quarks gain dynamical masses and become constituent quarks which 
are the valence quarks in potential model. The Goldstone boson pions are 
still massless. To account for this dynamical chiral symmetry breaking in our 
model, the coupling potential V must be large. This is why we introduce a 
constant term into the vector kernel. In this way, our model differs from the 
ordinary potential model (see ref. [1]) which had the constant term in the 
scalar interaction. In an ordinary potential model, this change will cause no 
difference (the potential U is not changed). But in this extended potential 
model, the coupling potential V can be made strong enough to give n mesons 
very low masses. 

Inserting the simple model wave functions (41) and (42) into Eq. (25) with 
the potentials given by (44) and (45), we obtain an eigen equation of a 2 x 2 
matrix 



A B 




X 




X 








= E 




-B -A 




Y 




Y 



(47) 



where 

^ 2 <^|) + ^-|-i(^>, <«) 
S = -f <<xi-<T 2 >. (49) 

The energy E can be easily obtained 

E = VA 2 - B 2 . (50) 



The wave function parameter (3 is determined by the minimum of the energy 
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.5 



0.5 











Fig. 1. Comparison of the integral / °° t 2 dtVt 2 + x 2 exp(— t 2 ) (solid line) in the av- 



erage quark kinetic energy with 1 + jx 2 (dashed line) used in Eq. (53). 
Since B is a constant, we get 



OA 



0. 



(52) 



We approximate the average quark kinetic energy by the following expression. 



(K) mc = (^k 2 + m 2 c }^ J-P 2 + m 2 c - 



(53) 



Numerically the deviation is very small and less than 5% (See Fig. 1). 

The current quark mass m can be obtained from the partial conservation of 
axial- vector current (PCAC). Consider the axial- vector current 



A» d (x) = u{x)^ lb d{x). 



(54) 



According to PCAC, its divergence is 

%Kd( x ) = ( m « + m d )u(x)i>y 5 d(x), 



(55) 
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where m„ and rrid are current masses of up and down quarks. Let us calculate 
the matrix element (0 | d^A^ d (x) \ n + ) for 7r + meson with momentum q n . One 
obtains 

(0 | dMx) I vr+, q„) = -ug(0 | A^(x) | 7T+, q w ). (56) 

/,From right hand side of PCAC relation (55), one also has 

(0 | dpA^Jx) | 7T+) = (m u + m d )(0 | u{x)i lb d{x) | tt+). (57) 

Both of the matrix elements can be easily calculated using the model wave 
functions for a static rr + meson. We obtain 

2m = m u + m d = *7 ^ Y? m n . (58) 



Before the numerical calculation, one can make some qualitative analysis of 
the properties of i\ and p mesons in the model. Let 



u = X + Y, (59) 

v — X — Y. (60) 

The matrix equation (47) can be written for u, v 

(A-B)v = Eu, (61) 

(A + B)u = Ev. (62) 

Also the normalization relation (30) becomes 

X 2 - Y 2 = uv = 2E. (63) 



First, let us consider the tt meson which is a Goldstone boson. We have 



f 16c\ 

2 A) g- J Vir = m^u*, (64) 

(^2A + ^ju 7T = m 7T v n , (65) 



where 
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If the current quark mass mo = 0, chiral symmetry is a strict QCD symmetry. 
From the Goldstone theorem, 7r meson should be massless. As an effective 
theory of QCD, the n meson mass in our model should also be zero. Assuming 
m c is the constituent quark mass in the limit of m = 0, we have 

2A)L c -^ = 0. (67) 



In the real world, current quarks have small nonzero masses. So does the ir 
meson. As the average mass of current quarks (u and d) slightly changes to 
m 7^ 0, the constituent quark mass will also change a little bit to m c + Am c 
with Am c <C m c . So we have 

4)L C+ Am c ~ A)L C =A) (68) 

and 

M^-M^^-A^^A-V (69) 



Eqs. (64) and (65) become 



Am c m c 

2—— v 7r = m 7r u 7T , (70) 

3AQU n = m n v 7T . (71) 
We obtain 



/ 6 Am c m c A 



With the normalization 2m n = u n v n , we obtain 



(72) 
(73) 



v v = yj6A^>. (75) 
/,From PCAC (eq. (58)), the mass of current quark is related to Am c as 



12 



Next we can calculate the weak decay constant of 7r meson to obtain 



(77) 



Eq. (73) means u n /v n ~ and Y n ps —X n , which shows the importance of the 
contribution to the it meson from the small component of y-channel . 



For the p meson, with the help of Eq. (67), the equations are 



2A v p = m p u p , 
A u p = m p v p . 



One simply obtains 



(78) 
(79) 



m p = V2A , 



u. 



V2. 



(80) 
(81) 



With the normalization 2m p = u p v p , 



u p = 2JA 



v p = J2A q . 

Next the electromagnetic decay constant of p meson can be calculated as 



(82) 
(83) 



Jo .__ 1 



TO, 



2Nc(jr 



(84) 



A relation among masses of n, p mesons and current quark mass m can be 
easily obtained 



m\ = Qm A = 3V2m m p . 



(85) 



^From the experimental p meson mass m p = 770 MeV, one can estimate that 
A ps m p j\[2 ~ 540 MeV. Then we can estimate the current quark mass to 
be toq ~ 6 MeV from the physical pion mass m n = 140 MeV. 
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Another relation about n and p decay constants can also be easily obtained 
from Eqs. (77) and (84), 



The KSRF relation [8,9] gives 

Uf P = ^| m p' ( 87 ) 



while the experiment value is m p /(f n f p ) = 1.62. In our approach Y p /X p = 
1/(3 + V2)= 0.17 from Eq. (81). This means the contribution to the rho 
meson from y-channel is rather small. If the small component is completely 
neglected, i.e., Y p = 0, we may obtain m p = |A and X p = y/'2m p . Then 



J p 



m p 



\ 




(88) 



thus 



Ufp = y| m P' ( 89 ) 



which is exactly the KSRF relation. The experimental value is between our 
result and KSRF relation. 

In numerical calculation we parametrize the running coupling constant as 



a s (Q 2 ) = $> fe exp(-Q 2 /4 7 2 ) 

k 

= 0.25 exp(-Q 2 ) + 0.15 exp(-Q 2 /10) + 0.20 exp(-Q 2 /1000) (90) 

as in Ref. [1] (with Q in GeV). Also the constituent quark mass m c is fixed to 
be m c = 220 MeV. The results are listed in Table 1. First (Set 1), we choose 
the confinement parameter b = 0.18GeV 2 according to lattice calculation. 
The p meson mass is somehow larger. Next (Set 2), we adjust the confinement 
parameter slightly to b = 0.15GeV 2 . Both of pion and rho meson masses agree 
with the experimental values well. The current and constituent quark masses 
are also reasonable in the ranges of theoretical and experimental estimations. 

The decay constants /V and f p deviate from the experimental data by a factor 
2. This is not very surprising. In the derivation of the effective Hamiltonian 
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Table 1 

Numerical results for it and p mesons. Experimetal data are taken from ref. [18]. 





b 


c 


P 


m 


m w 


m p 




(GeV 2 ) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


Exp. 










139 


770 


Set 1 


0.18 


220 


391 


6.5 


152 


843 


Set 2 


0.15 


OA A 

204 


364 


5.6 


136 


779 




Y-rr 


Yn 


U 
(MeV) 


f P 


Z 




Exp. 






93 


5.1 






Set 1 


-0.842 


0.169 


190(93) 


2.6(5.3) 


0.70 




Set 2 


-0.847 


0.169 


178(89) 


2.5(5.1) 


0.71 





(which we skip in this work), the current quark field \I/ is renormalized into 
constituent quark field. This means that the constituent quark may include the 
contributions from configurations other than the single quark, like qqq, qg etc 
as indicated through the deep inelastic scattering experiments. For example, 
we may assume the normalized constituent quark field ^/^ 

^jv = Zty + contributions from qqq, qg, .... (91) 

Here we have introduced the factor Z. Roughly speaking, \Z\ 2 is the proba- 
bility of finding the "pure" current quark field in a complicated constituent 
quark. However, from the electro-weak Lagrangian we know that what really 
participates in the electro-weak interactions are the current quarks. In other 
words, only the "pure" current quark will contribute to the pion weak decay 
constant and rho meson electromagnetic decay constant. In our calculation 
we need two real current quarks to annihilate. So there will be a difference 
of factor Z 2 . Decay constants f n and f p can then be fit to the experimental 
values if we set the additional normalization factor Z m 0.7 (In Table 1 the 
results are listed in parentheses). 



4 Summary 

The chiral symmetry is studied in an extended constituent quark potential 
model which includes the small component effect with RPA approximation. 
With some modification of quark interaction in potential model, 7r meson 
is still a Goldstone boson. The small mass of tc meson can be connected to 
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nonzero current quark mass from PCAC. All other mesons such as the vector p 
meson can also be studied in this model like in the ordinary constituent quark 
models. With the small component effect of y-channel contribution, the well- 
known KSRF relation of the ir and p meson decay constants is reestablished 
in the model. 

This extension to constituent quark potential model aims at improving the 
constituent quark model to comprise the small component, i.e., the y— channel 
quark excitations which are essential to rr (and K) mesons. To thoroughly 
investigate the effect of this extension, we need perform a full calculation of 
meson spectroscopy. 
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